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INTRODUCTION 


It  is  well-known  that  lower  ionospheric  layers  such  as  the  D- layer 
and  the  stratospheric  ozone  layer  exhibit  large  variations  of  density  (see 
Fig.  1),  which  are  reminiscent  of  gravity  wave  structures.  Further,  in 
some  thin  layer  structures,  such  as  noctilucent  clouds  (Fogle  and  Haurwitz, 
1966)  and  the  sporadic-E  layer  (Miller  and  Smith,  1975;  1976),  the  existence 
of  wave-like  structure  has  become  a singularly  dominant  feature.  Since  the 
presence  of  gravity  waves  cf  low  frequency  and  short  vertical  wave-length 
is  expected  in  the  atmosphere,  the  exhibition  of  wave  structure  in  atmo- 
spheric and  lower  ionospheric  layers  is  not  at  all  surprising.  What  is 
surprising,  however,  is  that  the  wave-like  density  variations  of  the  layer 
structures  are  generally  of  large  amplitude  whereas  the  density  variations 
in  the  background  atmosphere  are  usually  observed  in  the  middle  and 
upper  atmosphere  to  be  Less  than  or  of  the  order  of  10-15%  (see  Fig.  2a,  b). 
For  example,  the  wave-like  partial  pressure  variations  of  the  ozone  layer 
shown  in  Fig.  1 are  of  order  of  50  p»mb  out  of  the  layer  background  partial 
pressure  of  some  100  ^imb.  Density  variations  associated  with  such  wave- 
like structures  in  the  middle  atmosphere  have  been  observed  (Faire  et  al.  , 
1974).  Bounds  on  the  stratospheric  wave-like  density  variations  can  also  be 
estimated  from  observed  wave-like  horizontal  wind  components  in  the  strato- 
sphere (Webb,  1966)  by  means  of  perturbation  theory  (Hines,  I960).  The 
result  is  that  wave-like  density  variations  in  the  stratosphere  are  less  than 
10%  of  the  background  stratospheric  density,  in  agreement  with  direct  density 
observations.  Indeed,  one  may  contrast  wave-like  density  variations  of  the 
D and  E regions  of  the  ionosphere  with  that  of  the  neutral  atmosphere  in  the 
same  altitude  region  and  find  a similar  relationship. 
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Figure  1.  Stratospheric  ozone  layer  showing  wave  structure. 
(Komhyr  and  Gross,  1968). 
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Figure  Za.  Wave  structure  in  the  stratospheric  total  density. 
(Faireetal.,  1974). 
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Figure  2b.  Wave  structure  in  the  lower  thermospheric  density 
measured  by  a satellite -borne  accelerometer 
(Ching,  unpublished  data). 
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While  numerous  authors  have  from  time  to  time  qualitatively  espoused 
the  view  that  wave- like  density  variations  observed  in  the  atmospheric  and 
lower  ionospheric  constituent  layers  are  somehow  due  to  layer  response  to 
gravity  waves  in  the  background  atmosphere  rather  than  due  to  some  chemical 
or  mechanical  instability  in  the  layer  itself,  the  very  large  amplitudes  observed 
in  the  density  fluctuations  of  the  layer  structures  remain  a major  obstacle  to 
a quantitative  resolution  of  such  a view.  Recently,  the  unifying  feature  of  the 
gravity  wave  response  mechanism  was  further  strengthened  by  a consideration 
of  Dudis  and  Reber  (1976)  who  showed  that  density  variations  in  the  minor 
constituents  of  the  neutral  thermosphere  obey  relationships  consistent  with 
gravity  wave  perturbation  of  the  background  atmosphere,  although  the  thermo- 
spheric collisional  frequency  may  not  be  large  enough  to  support  the  approxi- 
mation that  the  constituents  are  strongly  coupled  to  the  wave  motion  by 
collisions  (Del  Genio  et  al. , 1977). 

In  this  paper  we  give  a calculation  of  the  linear  and  quasi-linear  re- 
sponse of  a constituent  layer  of  general  density  profile  to  a gravity  wave  per- 
turbation of  the  background  atmosphere.  It  is  assumed  that  the  motion  of  the 
layer  is  collisionally  coupled  to  the  wave  motion  of  the  background  neutral 
atmosphere;  therefore,  our  results  are  applicable  to  layers  below  200  km 
altitude  only.  Above  200  km  altitude,  the  layer  motion  depends  on  gradient  and 
diffusive  drifts;  however,  since  there  are  no  interesting  layer  structures, 
except  for  the  F- layer  of  the  ionosphere,  above  200  km,  our  assumption  is 
quite  adequate  for  the  purpose  of  considering  layer  responses.  We  shall  show 
that  the  amplitude  of  the  density  response  of  the  layer  is  a strong  function  of 
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the  gradient  of  the  layer  density  profile  relative  to  the  scale  height  of  the 
background  neutrail  atmosphere.  For  sharply  defined  layers  such  as  the 
ozone  layer  and  the  D and  E regions  of  the  ionosphere,  the  density  re- 
sponse of  the  layer  is  of  considerably  larger  amplitude  than  the  wave 
density  variation  in  the  background  neutral  atmosphere.  For  cases  where 
the  layer  response  can  no  longer  be  regarded  as  a small  perturbation, 
the  exact  quasi-linear  solution  shows  singular  behaviour  at  sufficiently 
sharp  layer  density  gradients,  thus  signifying  the  onset  of  turbulence 
forced  by  gravity  waves. 
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RESPONSE  TO  PERTURBATIONS. 


IL 

Consider  a constituent  layer  of  density  n(z)  embedded  in  the  neutral 
atmosphere  of  background  density  Nq(z)  which  is  stratified  with  scale  height 
Hq.  a linear  gravity  wave  perturbation  of  density  N,  horizontal  wind  com- 
ponent V and  vertical  wind  component  V is  assumed  to  exist  in  the  back- 

X z 

ground  atmosphere.  According  to  linear  gravity  wave  theory,  V and  V 

X z 

are  related  to  N/Nq  by 

= ( ikg/uu)  ( N/Nq  )[tD^YHQ-g(Y-l)]/[iUJ^+g(Y-l)K2.^ 

= ( N/Nq  ) [ - gtD(  1 - IyHqK^  ) ] / [ itu^  + g ( y - 1 ) ] (2) 

where  g is  the  gravitational  acceleration  and  Y = c /c  , the  ratio  of  heat 

p V 

capacities.  The  wave  density  perturbation  N/ Nq  is  given  in  terms  of 
amplitude  A and  phase  <{>  by 

N/N^  = A exp  ^ = A exp  ( iuut  - ikx  - iK  z ) (3) 

u z 

where  K = k + i/ZH^..  Equation  (3)  defines  the  wave  frequency  uu  , 

z z u 

horizontal  wave  number  k and  vertical  wave  number  k • The  gravity  wave 

z 

perturbation  forces  a response  in  the  constituent  layer  density  n and  the 
constituent  velocity  v . For  a neutral  constituent  layer,  the  constituent 
momentum  equation  relates  the  time  rate  of  wave  velocity  variation  to  wave 
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gradient  drifts  and  to  collisional  coupling  to  the  background  wave  velocity  V, 
Since  the  neutral  atmospheric  collisional  frequency  below  200  km  is  orders 
of  magnitude  larger  than  the  wave  variations,  the  constituent  momentum 
equation  reduces  to 


V = V , (4) 

which  is  the  usual  assumption  that  the  layer  moves  with  the  background  medium. 
For  an  ionospheric  layer  in  the  E and  FI  regions,  the  collisional  coupling  at 
frequency  between  ions  and  neutrals  and  the  gyration  about  the  magnetic 

field  with  frequency  = eB/m^^c  are  both  much  larger  than  the  gravity  wave 
frequency  and  drifts;  therefore  the  constituent  momentum  equation  is  approxi- 
mately 

X - 0 . (5) 

for  low  frequency  gravity  wave  perturbations  in  the  lower  ionosphere.  Inclusion 
of  gradient  drift  terms  in  (4)  and  (5)  can  be  easily  accommodated;  however, 
these  terms  will  detract  from  the  simplicity  of  our  analysis  without  adding  any 
significant  content. 

For  constituent  layers  whose  chemical  and/or  ionization  equilibrium 
lifetimes  are  long  compared  to  the  gravity  wave  period, such  as  the  ozone  layer 
and  the  sporadic- E layer,  the  response  of  layer  density  n is  governed  by  the 
continuity  equation 
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^ + V.  ( vn) 


0 


f 


(6) 


where  v is  related  to  the  gravity  wave  forcing  velocity  V by  either  (4)  or  (5), 
For  ordinary  ionospheric  layers  in  the  D and  E regions  at  night,  the  recom- 
bination lifetimes  are  comparable  to  gravity  wave  periods;  therefore,  the 
effects  of  recombination  must  be  taken  into  account  in  (6).  The  realization  of 
the  importance  of  recombination,  when  taken  together  with  the  results  of  the 
present  analysis,  is  crucial  in  the  consideration  of  nighttime  El- region  tur- 
bulence; however,  for  the  sake  of  simplicity  these  effects  will  be  dealt  with 
elsewhere.  Hence,  our  considerations  here  as  applied  to  layers  whose  re- 
combination lifetimes  are  shorter  than  or  comparable  to  gravity  wave  periods 
will  be  regarded  as  indicative  rather  than  quantitative  results.  For  the  re- 
mainder of  the  analysis,  results  based  on  (5)  will  be  given  without  proof, 
since  the  derivation  is  entirely  analogous  to  the  one  presented  in  terms  of  (4). 

Without  making  a small  perturbation  assumption  for  layer  response, 
the  constituent  layer  density  n can  be  written  as  the  sum  of  a layer  density 
profile  in  the  absence  of  a gravity  wave  in  the  background  atmosphere  nQ(z) 
and  a gravity  wave  density  response  n ( x,  z,  t ),  i.  e. 

n = Hq  (z)  + n ( X,  z,  t ) (7) 

In  our  calculation,  nQ  (z)  is  regarded  as  a general  layer  density  profile  generated 
by  some  production  mechanism.  Substitution  of  (4)  and  (7)  into  (6)  yields  the 
partial  differential  equation 
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(8) 


dn 

+i;.V)(S/no)+(^3#-  iKJ(l+S/no)V^.ik(l+S/n„)V^  = 0, 

which  relates  the  layer  density  response  (n/n^)  to  the  background  gravity  wave 
perturbation  ( N/Nq  ) via  (1)  auid  (2),  The  layer  density  profile  n^  manifests 
itself  in  (8)  in  terms  of  a gradient  length  function  L(z),  defined  as 

-1 


1 ^^0 
L(z)  = [—  ^ ] 

' ' n^  dz  •' 


(9) 


The  boundary  condition  defined  for  (8)  is  that  (n/n^)  should  vanish  for  vanishing 
gravity  wave  amplitude  A in  ( N/Nq)  = Ae^.  The  general  solution  to  (8)  is 
very  difficult  to  obtain  not  only  because  (8)  depends  on  "second  order"  terms  of 
the  form  {u/iIq)  ( N/Nq  ) but  also  because  L is  an  arbitrary  function  of  z. 
Thus,  the  exact  vertical  structure  of  the  response  requires  specification  of  a 
model  of  n^.  However,  for  purposes  of  interpretation  it  is  much  more  advan- 
tageous to  consider  the  local  approximation  for  which  L(z)  is  a slowly  varying 
function  since  the  background  wave  is  almost  never  a steady  wave  train.  In  this 
limit,  it  is  a simple  matter  to  solve  the  linear  response  version  of  (8).  A more 
exact  solution  of  (8)  can  be  sought  in  the  quasi-linear  limit  that  the  response 
(n/nQ),  similar  to(N/NQ),  depends  on  (x,  z,t)  through  the  phase  function 
0 = i uut  - ikx  - iK^z.  Inspection  of  (8)  shows  that  indeed  (n/nQ)  is  purely  a 
function  of  0 if  the  gradient  length  L(z)  is  large  compared  to  ( 8 0 / 8z  ) , 
i.  e.  k >>  1/L(z).  Strictly  speaking  the  quasi-linear  limit  is  invalid  at  the 
layer  peak  where  dnQ/dz  vanishes;  however,  separate  investigation  in  the 
neighborhood  of  the  layer  peak  shows  that  the  quasi-linear  solution  given  be- 
low approaches  the  proper  small  perturbation  response  even  when  Ij(z) 
approaches  infinity,  i.  e.  when  dnQ/dz  vanishes.  In  the  quasi-linear  limit, 

(8)  can  be  transformed  into  a non-linear  ordinary  differential  equation  in 
0 , whose  solution  is  given  by 


( n /n.  ) = ( a/  6)  ^ [ 1 + ( 6 / 8)  ( N/N.  ) ] 


(10) 
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where 


a=  - uu^  + ( g/Hg  ) ( 1 + Hq/L  ) - igK^  ( 1 + YHq/L  ) 
p = - uu  ^ + ig  ( Y - 1 ) K 

6=  - ( g/L  ) ( 1 - i Y K ) + 2 ou^  + 2ig  ( Y - 1 ) K 

U z z 

2 

For  low  frequency  gravity  waves,  i.  e.  ou  <<  g (Y-  1 )/(YHq)  , 

and  ou^k^  k^g  ( Y - 1 ) /(  Y ) , we  obtain 
z u 

[l+VHg/L]  , 

(n/no)~-YTr  ■ [i+YHp/jV-'l  t'L]  • »"  { 1 + ( WN„  ) [ 2 +VHo/(  V ■ 1)  L ] | (11 

In  the  small  perturbation  limit,  ( N/Nq  ) is  infinitesimal  and  (11)  approaches 
the  limit 

( n / nQ  ) - - ( 1 + YHq/L  ) ( N/Nq  ) , (12) 

which  is  exactly  the  result  obtained  from  (8)  under  the  linear  approximation. 
Further,  for  L = - , the  constituent  scale  height  of  a stratified  atmospheric 

constituent,  (12)  is  exactly  the  result  of  Dudis  and  Reber  (1976). 

For  the  ionospheric  case  (5),  the  situation  becomes  more  complex  since 
the  orientation  of  the  magnetic  field  is  involved.  If  we  consider  equatorial  E 
and  FI  regions  in  the  case  for  which  the  magnetic  field  is  assumed  to  be  per- 
pendicular to  the  X - z plane,  an  analogous  quasi- linear  solution  to  (6)  can 
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again  be  found,  although  the  validity  of  qua  si- linearity  depends  also  on 
k >>(  d n./dz  )/  Q.  and  k >>  ( d v.  /dz  )/  v.  « The  exact  quasi-linear 
solution  for  this  case  is  too  complicated  to  be  exhibited  here;  instead,  the 
linear  limit  analogous  to  (12)  will  be  given  as  a basis  for  interpretation  in 
the  next  section. 

(n./ng)  --  {l  + ivYHQk^+(YHQ/L)[l- v(l  + ak/v.^k^)]  } (N/Nq)  (13) 

where 


v = ( k^/k  ) [ ( ^in  ^ ^ 

Hi  the  derivation  of  (13)  and  (14),  the  low  frequency  gravity  wave  approximation 
has  been- assumed.  Note  that  in  the  D-region  >>  fi  and  v <<  1,  thus 

the  response  behaves  exactly  as  that  of  a neutral  layer.  In  the  F- region 
v.^  <<  and  (13)  approaches  the  limit  ( n^/ng  ) »-  ( N/Nq  ) /(  y - 1 ) , 
which  does  not  depend  on  the  layer  structure. 
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INTERPRETATION 


In  this  paper  we  shall  limit  our  interpretation  of  the  results  of  the 
previous  section  in  terms  of  their  application  of  layer  structures  in  the 
atmosphere  and  lower  ionosphere.  The  interpretation  of  the  linear  results 
for  the  particular  case  of  density  relations  between  stratified  neutral  thermo- 
spheric constituents  of  different  masses  has  been  undertaken  by  Dudis  and 
Reber  (1976). 

For  a general  layer  structure  such  as  the  ozone  layer  shown  in 
Fig.  1,  the  inverse  gradient  length  1/L(z)  is  positive  in  the  bottomside, 
vanishes  at  the  layer  peak  and  is  negative  in  the  topside.  According  to  (12), 
which  is  more  transparent  for  interpretation  purposes  then  (10)  or  (11),  one 
immediately  obtains  the  result  that  the  bottomside  response  is  necessarily  of 
larger  amplitude  than  the  topside  response,  whereas  at  the  layer  peak  the 
response  is  given  by  ( n/nQ)  =“  - ( N/Nq  ) / ( y-  1 )•  Such  a relationship 

is  generally  borne  out  in  an  examination  of  ozonesonde  observations  (Komhyr 
and  Sticksel,  1967;  Komhyr  and  Grass,  1968).  The  magnitude  of  response 
( n / nQ  ) relative  to  the  background  wave  ( N/Nq  ) depends  on  the  magnitude 
of  the  layer  profile  gradient  L(z)  relative  to  the  atmospheric  scale  height  Hq. 
For  a sharp  layer  density  profile,  | yH^/Li  [ > 1,  the  layer  density  re- 
sponse  n/nQ  can  be  considerbly  larger  than  the  background  gravity  wave 
density  variation  N/Nq.  For  example,  on  the  bottomside  of  the  ozone 
layer  shown  in  Fig.  1 the  "amplification”  factor  | ( 1 + y )/(  Y - 1 ) | 

is  approximately  5,  whereas  on  the  topside  this  factor  is  approximately  2. 
Further,  it  is  of  interest  to  note  that  if  | y Nq/L  | > 1 then  the  phase  of  the 
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response  ( h/hq  ) is  180°  out  of  phase  with  ( N/Nq)  at  the  bottomside  and  near 
the  layer  peak,  but  the  phase  difference  takes  a 180°  jump  just  above  the  layer 
peak  such  that  ( n / Uq  ) becomes  in  phase  with  ( N/Nq)  at  the  topside.  The 
apparent  180°  phase  jump  of  ozone  density  variations  just  above  the  layer 
peak,  shown  in  Fig.  1,  may  perhaps  be  interpreted  in  terms  of  this  relation- 
ship. In  the  bottomside  O and  E regions  of  the  ionosphere,  very  large 
wave-like  density  variations  have  been  observed  ( e.  g.  Pfister  and  Ulwick, 

1958;  Knight,  1972).  Again,  these  observations  are  qualitatively  consistent 
with  the  relationship  between  the  "amplification"  factor  of  (12)  or  (13)  and 
the  variation  of  the  layer  density  gradient  L(z),  although  our  results  must  be 
modified  by  inclusion  of  recombination  effects. 

The  sporadic- E layer  of  metallic  ions  is  of  particular  interest  since 
the  layer  is  sharply  defined  and  the  recombination  lifetime  is  long  compared 
to  gravity  wave  periods.  Recently,  incoherent  backscatter  experiments  have 
shown  that  the  layer  is  strongly  turbulent  and  exhibits  large  variations  of 
density  (Miller  and  Smith,  1975;  1976).  Indeed,  at  times  definite  gravity  wave 
effects  can  be  identified  (Miller  and  Smith,  1976),  We  wish  to  point  out  here 
that  the  "amplification"  factor  for  the  sporadic- E layer  can  be  quite  large. 

At  the  vicinity  of  100  km  altitude  where  sporadic-E  occurs  v fi.  ~ 30; 
therefore  the  ionospheric  layer  relation  (13)  collapses  to  that  of  the  neutral 
layer  relations  (11)  and  (12).  Since  L for  the  sporadic-E  layer  is  of  order 
of  1 km  and  Hq  is  of  order  of  8 km,  yH^/Li  is  sufficiently  large  that  the 
linear  relation  (12)  becomes  invalid  and  our  interpretation  must  rely  on  the 
quasi-linear  relation  (11),  For  a background  gravity  wave  with  | N/Nq  | - 10%, 
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we  obtain  the  sporadic-E  layer  response  cunplitude  | n/  nQ  | 1.4, 

which  certainly  signifies  the  formation  of  very  dense  knots  of  sporadic-E 
layer  density  observed  by  Miller  and  Smith  (1976),  Indeed,  for  gravity 
waves  with  k — 1/L(z),  even  the  quasi-linear  solution  becomes  invalid 
and  the  response  ( n / n^  ) is  dominated  by  a resonance  between  the  layer 
structure  and  the  wave  structure  which  would  result  in  turbulence.  Although 
the  exact  non-linear  solution  for  such  a case  is  too  difficult  to  obtain,  that 
the  layer  response  presages  turbulence  in  such  a case  can  be  argued  from  a 
consideration  of  the  relative  phase  between  ( n / n^  ) and  ( N/Nq)  in  (11). 
When  I I ^ ^ (11)»  the  complex  argument  of  the  logarithm  in 

(11)  may  approach  its  negative  real  axis  branch  cut;  as  a result,  the  phase 
of  ( n / Uq  ) may  take  on  an  infinite  number  of  branch  values  separated  by 
multiples  of  n.  The  major  characterization  of  turbulence  is  that  the  phases 
of  its  density  variations  are  not  well  correlated.  Based  on  these  consider- 
ations, we  would  like  to  suggest  that  gravity  wave  forcing  is  the  major 
mechanism  in  the  generation  of  density  irregularities  in  sharply  defined  con- 
stituent layers  of  the  atmosphere  and  ionosphere.  Further,  sharp  layer 
gradients  are  necessarily  "turbulent". 

The  generation  of  turbulence  is  of  great  interest  in  the  nighttime 
equatorial  ionosphere  in  the  200  km  region.  It  has  been  suggested  that  the 
nighttime  equatorial  F- region  is  susceptible  to  Rayleigh- Taylor  instability 
(e.  g.  Woodman  and  LaHbz,  1976;  Hudson  and  Kennel,  1975);  however,  the 
Rayleigh- Taylor  growth  rate  in  the  200  km  region,  g/ ( 2 v ) * 10“^ 
sec“^,  is  sufficiently  small  that  gravity  wave  response  effects  at  wave 
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period  short  compared  to  Rayleigh- Taylor  growth  time  may  be  more  important. 
We  wish  to  point  out  here  that  the  gravity  wave  response  at  the  bottom  side  of 
the  layer  is  subject  to  Rayleigh- Taylor  instability  also.  An  analysis  of  the 
gravity  wave  response  of  an  ionospheric  layer  structure  in  the  limit 
shows  that  internal  wave  modes  are  also  susceptible  to  Rayleigh- Taylor 
growth.  Since  internal  wave  modes,  unlike  the  surface  waves  of  the  pure 
Rayleigh- Taylor  instability,  can  propagate  throughout  the  layer  structure,  we 
would  expect  that  gravity  wave  response  at  the  steepened  bottomside  of  the 
nighttime  equatorial  ionospheric  layer  structure  will  play  a major  role  in  the 
"bottomside  turbulence"  phenomenon.  Detailed  analysis  of  such  a situation 
will  be  given  elsewhere. 

As  a final  example,  we  wish  to  consider  the  atomic  oxygen  layer  which 
maximizes  slightly  above  the  mesopause  at  altitudes  between  90  and  100  km. 
Although  very  few  direct  measurements  of  the  atomic  oxygen  distribution  in 
this  region  of  the  atmosphere  are  available,  some  examples  of  wave  structure 
may  be  seen  in  the  rocket-bome  mass  spectrometer  data  presented  by  Philbrick 
et  al.  (1973,  1974).  In  these  cases  the  wave  amplitudes  were  relatively  modest 
( '^lO  - 15%  ) and  no  wave  features  were  apparent  in  the  simultaneously  measured 
background  gas  profile.  This  behavior  is  consistent  with  our  theory  since  the 
amplification  factor  is  predicted  by  (12)  to  be  about  2.5  at  the  layer  peak  and 
larger  below  the  peak.  Further  evidence  of  irregular  structure  in  the  atomic 
oxygen  distribution  near  the  layer  peak,  which  would  appear  as  wave-like 
os:cillations  in  a vertical  profile,  may  be  found  in  latitude- altitude  contours 
derived  by  Wasser  and  Donahue  (1978)  from  5577- A airglow  measurements.  In 
these  remotely- sensed  data  under-dense  "patches"  are  seen  particularly,  though 
not  exclusively,  on  the  bottomside  of  the  layer.  The  latitudinal  extent  of  the 
patches  is  a few  degrees,  i.  e.  of  the  order  of  a few  hundred  kilometers.  It 
would  be  interesting  to  determine  in  the  future  whether  wave  perturbations  of  the 
stratospheric  ozone  layer  may  similarly  be  "patchy"  in  their  horizontal  extent. 
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CONCLUSION 


By  comparing  observations  with  the  solution  of  the  linear  and  quasi- 
linear  problems  of  the  density  reponse  of  a neutral  or  ionospheric  layer 
structure  of  gravity  wave  perturbations  in  the  background  atmosphere  below 
200  km  altitude,  we  conclude  that  large  wave- like  and  turbulent  density 
variations  in  the  ozone  layer  and  the  sporadic- £ can  be  interpreted 
as  layer  density  response  to  gravity  waves.  Limited  observations  of  wave 
structure  in  the  atomic  oxygen  profile  near  its  maximum  also  support  this 
interpretation.  The  layer  density  response  is  generally  of  larger  amplitude 
at  the  bottomside  than  at  the  topside,  although  the  amplitude  of  both  depends 
on  the  relative  sharpness  of  layer  density  gradient  in  comparison  to  the 
gradient  of  background  atmosphere  stratification.  It  is  suggested  that 
gravity  wave  response  plays  a major  role  in  the  generation  of  turbulence 
in  the  nighttime  equatorial  ionosphere  at  the  200  km  region. 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  Investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reenti^  physics,  chemical  kinetits,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high- power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena.  Including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro -optics;  conrununi cation  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far -infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  oi  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry:  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  air  glow;  magneto  spheric  physics,  cosr^c  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphero,  and 
magnetosphere:  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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